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Two-dimensional (2D) materials are not expected to be 
metals at low temperature due to electron localization (T). 
Consistent with this, pioneering studies on thin films re¬ 
ported only superconducting and insulating ground states, 
with a direct transition between the two as a function 
of disorder or magnetic field |2-6). However, more re¬ 
cent works have revealed the presence of an intermedi¬ 
ate quantum metallic state occupying a substantial region 
of the phase diagram |7-10| whose nature is intensely 
debated ifTTUTTl . Here, we observe such a state in the 
disorder-free limit of a crystalline 2D superconductor, pro¬ 
duced by mechanical co-lamination of NbSe 2 in inert at¬ 
mosphere. Under a small perpendicular magnetic field, we 
induce a transition from superconductor to the quantum 
metal. We find a unique power law scaling with field in 
this phase, which is consistent with the Bose metal model 
where metallic behavior arises from strong phase fluctua¬ 
tions caused by the magnetic field lHIiflD . 

Global superconductivity emerges in a sample when con¬ 
duction electrons form Cooper pairs and condense into a 
macroscopic, phase-coherent quantum state. In two dimen¬ 
sions, the phase coherence can be disrupted even at zero 
temperature by increasing disorder, either by degrading crys¬ 
tal quality or applying magnetic fields to create vortices (21. 
Granular or amorphous superconducting thin films, for which 
disorder levels can be controlled during growth, have thus 
provided an established platform for the study of quantum 
phase transitions in 2D superconductors. Within the conven¬ 
tional theoretical framework, increasing disorder or magnetic 
field perpendicular to a strongly disordered film at T = 0 in¬ 
duces a direct transition to an insulating state as the normal 
state sheet resistance approaches the pair quantum resistance 
h/( 2e) 2 = 6.4 kO (210). As film quality has improved over 
time, however, an intervening metallic phase with resistance 
much lower than the normal state resistance has been observed 
in several systems with generally less disorder (TiflOl. Its ori¬ 
gin is not well understood, and the various theoretical treat¬ 
ments can be generally divided between purely bosonic-based 
models, in which Cooper pairing persists in the metallic phase 
but phase coherence is lost 03HHI, and models that also in¬ 
corporate other fermionic degrees of freedom lfl5ifT7l . 


Recently, mechanical exfoliation has emerged as a tech¬ 
nique to produce ultra-clean, crystalline 2D materials, with 
graphene being a well-known example d). Like amorphous 
films, the thickness of these samples can be easily controlled 
down to the level of individual atomic layers. In contrast to 
amorphous films, a 2D superconductor exfoliated from a lay- 



Figure 1 . Environmentally controlled device fabrication, a) 

Schematic of heterostructure assembly process. Boron nitride 
(BN)/graphite (G) on a polymer stamp (PDMS) is used to electri¬ 
cally contact and encapsulate NbSe2 in inert atmosphere. The het¬ 
erostructure is lithographically patterned and the edge of graphite is 
metallized with Cr/Pd/Au. b) Optical images of heterostructure be¬ 
fore (left) and after device fabrication (right). In the (false-colored) 
left panel, the bilayer NbSe2 is outlined in a dashed green line and 
the overlap between the graphite and bilayer NbSe2 is shaded black. 
Scale bar is 5 pm in both images. 
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ered, single crystal, can exist in the regime of minimal disor¬ 
der, allowing for new insight into the nature of the vortex state 
in two dimensions. In this work, we realize such a system 
using a clean bilayer of NbSe 2 , a well-known type-II super¬ 
conductor with T c ~ 7.2 K in bulk form [19’ 201. Unique to 
this sample, the normal state sheet resistance is two orders of 
magnitude below h/ (2e) 2 and insulating behavior is never ob¬ 
served. Instead, the intermediate metallic phase emerges in an 
exceptionally large region of the magnetic field-temperature 
phase diagram. Unlike the typical exponential behavior asso¬ 
ciated with the quantum tunneling of fermionic quasiparticles 
mm , we observe a new power-law scaling as a function of 
field at low temperature that is consistent with the Bose metal 
scenario of the metallic phase [11 -141. 

Initial studies on exfoliated NbSe 2 flakes do not observe 
a superconducting transition to a zero-resistance state in the 
atomically thin limit 121II22II . Recently, however, it has been 
shown the surfaces of metallic materials may oxidize, alter¬ 
ing the electronic properties of thin samples [ 23 J. Exfoliation 
and encapsulation by a protective layer in inert atmosphere is 
thus crucial for preserving the intrinsic properties of the 2D 
material (231124) . We achieve this using a mechanical transfer 
setup installed inside a nitrogen-filled glove box (see Meth¬ 
ods). In short, within the glove box, an exfoliated NbSe 2 flake 
is electrically contacted by graphite (G), and the entire device 
is protected by an insulating layer of hexagonal boron nitride 
(BN). The graphite leads are then contacted using an edge- 
metallization technique in the ambient environment [25], 261. 
A schematic depicting the assembly and fabrication process is 
shown in Figure la and optical images of the heterostructure 
are shown in Figure lb before (left) and after device fabrica¬ 
tion. 

In the main panel of Figure 2a, we show four-terminal sheet 
resistance as a function of temperature for a particular NbSe 2 
device prepared using the method described above. The 
NbSe 2 thickness is 1.5 nm as determined by an atomic force 
microscope (AFM) outside the glove box after BN/graphite 
transfer, suggesting it consists of only two atomic layers. The 
resistance in the normal state is Rjy = 75 U. This corresponds 
to a residual resistivity that is 10 times larger than that of 
bulk crystals lf20l . yet the sheet resistance is still an order of 
magnitude less than that of the most conductive amorphous 
superconducting films (Toll . We observe a clear supercon¬ 
ducting transition to a zero-resistance state measured to the 
limit of our instrument resolution. The critical temperature 
is T c = 5.26 K, as defined by where the resistance is 90% 
of the normal state value, which is slightly reduced from the 
transition temperature of bulk samples (7.2K) (191120) . 

In order to further characterize the quality and dimension¬ 
ality of our sample, we have measured the temperature- de¬ 
pendent critical fields, defined by R(H c2: T) = 0.9 Rn, for 
field configuration both perpendicular and parallel to the lay¬ 
ers, and the results are plotted in the inset of Figure 2a. Close 
to T c , we expect = ($o/27r£o)(l — T/T c ), where <I>o = 
h/2e is the flux quantum and £o is the in-plane coherence 
length at zero temperature. A linear fit shown by the black line 




Figure 2. Characterization of bilayer NbSe2 device, a) Sheet re¬ 
sistance with temperature shows superconducting transition at T c = 
5.26 K. Temperature-dependent critical magnetic fields parallel and 
perpendicular to the layers is shown in the inset. Black line is lin¬ 
ear fit to H ^2 oc 1 — T/T c at high temperatures. Red line is fit to 
h [2 oc yjl — T/T c , the Tinkham formula for 2D samples [27]|. b) 
Voltage-current behavior at different temperatures. Inset shows ex¬ 
ponent a vs. T extracted from power law fitting V ~ I a near the 
normal state transition, a - 3 at the BKT temperature 5.01 K. 


yields £o = 8.9 nm, similar to the bulk value (19] l20) . From 
the normal state resistance and carrier concentration as deter¬ 
mined by Hall measurements, we estimate the electron mean 
free path to be £ = 17 nm, smaller than in bulk crystals ['20], 
but nearly twice the coherence length, confirming that our de¬ 
vice is in the pure superconductor regime (27). In contrast, 
evaporated films are generally characterized as “dirty” super¬ 
conductors with £ 0 > L The critical parallel field does not ex¬ 
hibit linear temperature dependence expected for anisotropic 
3D superconductors f27) . Instead, a 2D superconductor with 
thickness d < £ obeys ffjj 2 = (vT2<f> 0 /27r£o d)^/l — T/T c 
127) . The red curve in the inset of Figure 2a shows a best fit to 
this expression, from which we extract d = 3.4 nm, which is 
slightly larger than the thickness as determined by AFM (1.5 
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Figure 3. Magnetic field tuned phase transitions in 2D NbSe 2 . a) 2D color plot of sheet resistance vs. temperature and perpendicular 
magnetic field, b) Arrhenius plot of resistance for several magnetic fields shows thermally activated regime (black lines) and saturation at low 
temperatures (colored lines), c) Energy barrier vs. magnetic field extracted from linear fit to activated region. Solid red line is empirical fit to 
formula in inset. 


nm). However, this fitting has been previously found to over¬ 
estimate the true sample thickness in sufficiently thin systems 

EH. 

The superconducting phase transition in a 2D material with 
d < £o is understood to be of the Berezinskii-Kosterlitz- 
Thouless (BKT) type (29). In this scenario, the low- 
temperature, zero-resistance phase consists of bound vortex- 
antivortex pairs created by thermal fluctuations. Upon heat¬ 
ing, the pairs dissociate and may move, inducing dissipation. 
The BKT temperature defines the vortex unbinding transition 
and can be determined using current-voltage measurements 
as a function of temperature T, as shown in the main panel 
of Figure 2b. Current excites free-moving vortices, causing a 
nonlinear voltage dependence: V I a ( T ) m At Tbkt, a 2D 
superconductor obeys the universal scaling relation, V / 3 
(solid line in Fig. 2b) (29ll30ll . In the inset, we plot a vs. T, 
as determined by the slope of the different V — I traces in 
log-log scale. An example guide-to-eye fit is marked by the 
dashed line in the main panel. We determine Tbkt = 5.01 K 
from where a = 3 interpolates, only slightly less than T c as 
defined above. This is consistent with the behavior of systems 
with normal state resistance much less than h/{ 2e) 2 , where 
Tbkt is expected to be very close to the mean-field transition 
temperature ED 

The measurements performed above confirm that our de¬ 
vice exhibits the characteristics of a true 2D superconductor. 
The extracted material parameters together with the low nor¬ 
mal state sheet resistance further places the high sample qual¬ 
ity within a previously unexplored regime. We next turn to 
the dependence of resistance on perpendicular magnetic field 
as the effect of vortices can now be cleanly separated from the 
low static disorder present in the sample. Shown in Figure 3a 
is a 2D color map of the four-terminal sheet resistance in the 
same device as a function of both temperature and magnetic 
field applied perpendicular to the layers. For clarity, tempera¬ 


ture traces for different field levels are shown in an Arrhenius 
plot in Figure 3b. For fields larger than ~ 3 T, the sample 
is in the normal state for all temperatures. Previous works on 
strongly disordered films observed a field-tuned transition to 
an insulating state of bosons EEl; however, one might not 
expect insulating behavior for finite samples of a highly con¬ 
ductive 2D system in the disorder-free limit ED. We have 
applied perpendicular fields as large as 14.5 T and still see 
metallic behavior in the temperature dependence (see Supple¬ 
mentary Information, Figure SI main panel). 

As one lowers the field to just below H c 2 , a resistance drop 
is observed upon cooling from the normal state. In this region 
of the H — T phase diagram, the device exhibits activated be¬ 
havior, as can been seen in the linear slope in the Arrhenius 
plot (black lines in Fig. 3b). Classically, dissipation in a su¬ 
perconductor in which resistance is less than the normal state 
value can be attributed to the motion of individual vortices 
(flux creep or flow) (271 . In a clean 2D system, we expect the 
dominant energy barrier to flux motion to be that of vortex- 
antivortex dissociation I f32l . Flux resistance then becomes 
thermally activated when the temperature is comparable to 
the barrier energy. We have determined the activation energy 
from the linear portion of the Arrhenius plot in Figure 3b for 
different magnetic fields, and the result is plotted in Figure 3c. 
The functional form is expected to be U(H) = Uo ln(iT 0 /iT), 
where Uo = ^d/ (2567t 3 A 2 ), the vortex-antivortex binding 
energy, and Ho H c 2 El El. A fit to this form yields Uo 
= 27.5 K and H 0 = 2 T for our device. Assuming the mag¬ 
netic penetration depth A is similar to the bulk value (230 nm) 
E3, we estimate Uo ~ 14 K, on the order of the experimental 
value. 

At lower temperatures, the resistance saturates to a level 
dependent on magnetic field (see colored lines in Fig. 3b), 
hallmark of the quantum metallic state that is the main subject 
of this report. Qualitatively similar behavior has also been 
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Figure 4. Emergence of the quantum metal, a) Magnetoresistance above the superconducting transition for different temperatures. The 
data scale to a power-law R ~ (H — H c0 ) a ^ and collapse onto a single curve in the quantum metallic phase below IK. a vs. T is plotted 
in inset, b) Full H — T phase diagram of bilayer NbSe 2 device. The red circles are the locus of points where R(H,T) = 0.9 Rn- The 
purple squares, dividing the quantum metal from the thermally-assisted flux flow (TAFF) regime, show the transition from activated behavior 
R ~ exp(A(i7)/T) to temperature-independent R = R(H), i.e. the intersection of the black and colored lines in Fig. 3(b). The blue 
triangles denote the boundary of the superconducting phase H c q(T), the field at which resistance increases above the noise floor value. This 
criterion is determined by when hysteresis vanishes in V — I measurements (see Fig. S4). 


observed in amorphous MoGe and Ta films with much larger 
normal state sheet resistances 130 . This effect cannot be un¬ 
derstood within a classical framework, in which we expect to 
recover superconductivity for T U. In the past, various 
theories have been advanced to explain the origin of dissipa¬ 
tion in this temperature regime for disordered films EHED- 
One can distinguish between the different models based on 
the magnetic field dependence of the saturated resistance at 
low temperatures. 

Shimshoni et al. consider a disordered supercon¬ 
ducting film as a percolating network of superconduct¬ 
ing islands within an insulating matrix fl5l . Field- 
induced vortices that tunnel across thin superconducting 
constrictions give rise to resistance when coupled to a 
fermionic bath. The resistance depends on field as R = 
h/(4e 2 )exp[CTr/2((h/e 2 )/R N )((H - H c2 )/H c2 )\, where 
Rn is the normal state resistance and C is a dimensionless 
constant of order unity. This expression finds good agreement 
with the measurements of Ephron et al. on amorphous MoGe 
films, in which Rn ^ 1 kD ta, but evidently provides a poor 
fit to our data (see Supplementary Information, Figure S2). 
Furthermore, since our device is over an order of magnitude 
more conductive, the closest fitting requires C = 0.14, an un- 
physically small value. 

Galitski et al. consider a “vortex metal” phase where field- 
induced vortices interact with electrically neutral spinons 
ira . The theory accounts for the large peak in low- 


temperature magnetoresistance observed in disordered InOx 
films before saturation to Rn at higher fields Col. In our de¬ 
vice, however, we observe a monotonic dependence of resis¬ 
tance on field and see no peak structure in magnetoresistance 
up to 7 T (see Supplementary Information, Figure SI inset). 

Das and Doniach, as well as Dalidovich and Phillips, report 
that a 2D system of interacting bosons may form a gapless, 
nonsuperfluid state in the limit of zero temperature, a phase 
which they term a “Bose metal” CDfll. They argue that 
the uncondensed Cooper pairs and vortices are responsible for 
the small resistance observed at small finite field C2 A mag¬ 
netic field introduces gauge fluctuations, which disrupts phase 
coherence and causes dissipation, a quantum analog of that 
caused by thermal fluctuations in the BKT transition. In the 
Bose metal model, resistance on the metallic side of the field- 
tuned transition can be described by R (H-H c o) 2p , where 
H c o is the critical field of the superconductor to Bose metal 
transition and v is the exponent of the superfluid correlation 
length, which diverges across the boundary as (H — H c0 )~ u . 

In the main panel of Figure 4a, we show a log-log plot of R 
vs H — H c o taken at several different temperatures. H c o is a 
small temperature dependent value which we determined us¬ 
ing a method that shall be described below. The linear scaling 
observed here suggests a power-law dependence on field. We 
have fit the data to the expression R ~ (H — H c o) a and the ex¬ 
tracted exponent a is plotted in the inset as a function of tem¬ 
perature. At high temperatures, but below T c , resistance in- 
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creases roughly linearly with field, as expected for unhindered 
flux flow: R ~ R n H/H c2 (a = 1) f27l . The gray line shows 
linear scaling as a guide-to-eye. As temperature is lowered, 
the field dependence becomes increasingly nonlinear and col¬ 
lapses onto a single curve below 1 K with a ~ 3. The red line 
is an empirical fit given by R[Q] = 5.44(iT[T] — 0.15) 3,21 , 
which shows excellent agreement with the data. This yields 
a critical exponent of v = 1.61. While previous measure¬ 
ments on MoGe films at the lowest accessible temperatures 
observe mostly an exponential dependence of resistance with 
field (7; 331, which is consistent with quantum tunneling of 
fermionized vortices ca, at very small fields for the MoGe 
film, however, the scaling obeys a power-law with unity ex¬ 
ponent [33]. The power-law scaling we observe in our sample 
consistently over the entire field range thus suggests that the 
fermionic tunneling undergoes a crossover to Bose metal be¬ 
havior in the limit of vanishing static disorder. We also found 
similar field scaling in the quantum metallic state for addi¬ 
tional ultrathin devices, whereas this phase is distinctly absent 
in the bulk crystal (see Supplementary Information, Figure 
S3). 

The critical field of the transition out of the true zero- 
resistance state H c o is difficult to determine directly from the 
linear resistance given the limited accuracy of our instruments. 
Recently, however, experiments by Qin et al. and Li et al. on 
disordered Ta films showed that H c 0 can be determined indi¬ 
rectly using current-voltage measurements 0 .34]. Hysteresis 
is observed in the V — I characteristics on the superconducting 
side, similar to that seen in underdamped Josephson junctions 
E3, which disappear at the onset of the metallic phase. We 
have performed V — I measurements on the same bilayer de¬ 
vice at T - 0.5 K for increasing magnetic fields (see Supple¬ 
mentary Information, Figure S4 main panel). Clear hysteresis 
is seen between current sweep up and down for low fields. In 
the inset, we plot the current hysteresis A I for where the volt¬ 
age jump occurs as a function of magnetic field. A I vanishes 
close to H c o = 0.175 T. This then allows us to identify the 
true superconducting phase as that for H < H c o . We have re¬ 
peated this measurement for several different temperatures in 
order to determine critical field as a function of temperature. 

Figure 4b shows a full H — T phase diagram for our de¬ 
vice. The blue triangles mark H c0 (T) and the red circles mark 
77^2 (T). The boundary between quantum metal and flux flow 
(purple squares) is defined by the intersection of the fits to 
the activated resistance and saturated levels in Figure 3b. The 
normal phase extends up to at least 14.5 T without the appear¬ 
ance of an insulating phase, in contrast to previous works on 
strongly disordered films. The high sample quality is made 
possible by our facile device assembly technique in inert at¬ 
mosphere, which demonstrates a new route for the production 
of 2D superconductors in the ultraclean limit. 

Methods. 

(1) Crystal synthesis. Polycrystalline NbSe 2 was made by 
heating stoichiometric amount of Nb powder (99.5%) and Se 
shots (99.999%) in evacuated silica ampoules. Single crystals 
of NbSe 2 were grown by using a vapor transfer method. 400 


mg of NbSe 2 powder and 80 mg of I 2 were sealed in a 23cm- 
long silica ampoule with a 1.13 cm 2 inner cross section. The 
charge was put in the hot zone of 850°C and the sink is in the 
cold zone with 750° C. After one week, all the poly crystals 
became single crystals, while most of the thick plates of single 
crystals were found in the hot zone. 

(2) Device assembly and fabrication. We have exfoliated 
ultrathin NbSe 2 flakes in a nitrogen-filled glove box contain¬ 
ing less than 2 ppm oxygen. Separately, we prepared thin 
hexagonal boron nitride (hBN) on a polydimethylsiloxane 
(PDMS) stamp covered with polypropylene carbonate (PPC), 
which we use to “pick up” two closely spaced graphite flakes 
(separation ^ 2/im). The hBN/graphite stack is then used 
to cover the NbSe 2 flake inside the glove box. While graphite 
makes electrical contact to NbSe 2 , hBN provides an insulating 
oxidation barrier. Subsequent lithography may then be used to 
define a four-terminal device, in which the gapped region be¬ 
tween the graphite leads form the channel. The graphite is 
then electrically contacted using an edge-metallization tech¬ 
nique 125, 26|. 
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Figure SI. Behavior at higher magnetic fields. Sheet resistance with temperature for several different perpendicular fields. No insulating 
behavior is observed up to 14.5T. Inset shows magnetoresistance for several different temperatures. No peak is observed up to 7T. 



Figure S2. Vortex tunneling scenario. Magnetoresistance at 0.25K fit to the formula in Shimshoni et al. describing quantum tunneling of 
vortices with several different values for constant C GU. Best fitting (C = 0.14) deviates significantly at higher fields. 
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Figure S3. Thickness dependence of metallic behavior, a) Arrhenius plot of normalized sheet resistance with different magnetic field levels 
for 2.5 nm (left) and 5 nm (right) thickness NbSe 2 devices. Both show metallic behavior at low temperature and intermediate fields, b) 
Same plot for bulk crystal shows absence of the quantum metallic state. Transition temperature is reduced for increasing magnetic field, c) 
Magnetoresistance below the superconducting transition for different thin flake devices at 0.25K scaled to a power-law R ~ (H — H c o) a . a 
vs. thickness is plotted in inset. 
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Figure S4. Field-tuned voltage-current characteristics. V — I traces at 0.5 K for several perpendicular field levels. Inset shows hysteresis 
between current sweep up and sweep down as a function of field. Hysteresis is observed in the superconducting phase, H < H c o =0.175 T, 
but disappears for the metallic phase, H > H c q. 





















